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Abstract 

We study the potential of a static quark anti-quark pair in the range 0.05 fm < 
r < 0.8fm, employing a sequence of lattices up to 64^. Lattice artifacts in potential 
and force are investigated theoretically as well as numerically and continuum 
quantities are obtained by extrapolation of the results at finite lattice spacing. 
Consistency of the numerical results with the form of scaling violations predicted 
by an analysis a la Symanzik is found. The scale ro/a is determined for the Wilson 
action up to /9 = 6.92. 
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1 Introduction 



The A-parameter of quenched QCD has been computed [H] by means of an iterative 
finite size scahng method which allows to connect long distance non-perturbative 
quantities with short distance renormalized couplings in a controlled manner p[. 
The result is the A-parameter in units of the low energy scale rg. The latter is 
equivalent to the force between static quarks at intermediate distances r = tq ~ 
0.5fm 0. 

Given this information, the perturbative expansion for the static potential 
becomes a parameter free prediction for the force between static quarks at 
short distances. 

On the other hand, the potential itself was the very first quantity to be com- 
puted by Monte Carlo simulation of the pure gauge theory [§] and it is therefore 
very interesting to compare perturbation theory with the non-perturbative force. 
Since the original work of Creutz in the SU(2) gauge theory, many computations 
have been performed in SU(3) but the short distance region has not yet 

been investigated convincingly. The reason is simple: by "the potential" we do 
of course mean the potential in infinite volume; in practice we know that for the 
force at distances up to r 0.5 fm, the deviation from infinite volume are small 



on an torus with L = 1.5 fm |jT6|. With such a value of L, an investigation of 
distances of r 0.05 fm with control over 0(a/r) lattice artifacts requires very 
large lattices. 

In the A^f = theory (for our purposes this means in the pure Yang Mills 
theory) large lattices may nowadays be simulated. Still it would be very costly to 
perform a computation of the force in the full range say 0.05 fm < r < 0.8fm, since 
large distances on fine lattices require also large time extents of the Wilson loops 
and - with presently known techniques - compute intensive smearing methods. 

To investigate this range we therefore separately consider two regions, 0.05 fm < 
r < 0.3fm and 0.2 fm < r < 0.8fm, where the first one needs very large lattices and 
small lattice spacings and the second one had been simulated before on coarser 



lattices |jT6[. The region of overlap serves for calibration. Our reference scale, tq, 



is computed in the second range and related to Vc ~ 0.26fm, which is accessible 
in both. 

Apart from the wish to test renormalized perturbation theory, an additional 
motivation for the computation presented here is that we can enlarge the region 
of bare couplings (3 for which rg/a is known. It is now extended up to /3 = 6.92. 
This opens up the possibility to investigate scaling violations over a larger range 
of lattice spacings, which is of interest in view of the unexpected results found in 



the 2-d sigma model There indications for scaling violations differing from 



1 



their form expected from an analysis a la Symanzik have been found. In our 
paper we shall first apply Symanzik's theory of lattice spacing effects to the static 
potential, closing a gap in the literature. We then discuss our numerical results 
for a-effects. A detailed comparison of the non-perturbative force and potential 
with perturbation theory will be presented in a separate publication. 

2 The scales ro, 

2.1 Definition 

In a pure gauge theory we have to specify one physical quantity in order to renor- 
malize the theory. While a renormalized coupling in a non-perturbative scheme 
may be used, it is more convenient to choose a dimensionful long distance observ- 
able, which introduces the scale into the theory. Predictions of all other dimen- 
sionful quantities are then expressed in units of this scale. They are well defined 
and may be extrapolated to the continuum limit from results at finite resolution. 
It is clear that the scale should be chosen with care since it influences the preci- 
sion of many predictions. The length scale ro, defined in terms of the force F{r) 
between static quarks by the implicit equation 

r^F(r)U(,)=c, ro = r(1.65), (2.1) 

has turned out to be a good choice: it may be computed with good statistical and 
systematic precision. In QCD, tq has a value of about 0.5 fm 0. 

As discussed in the introduction, it is convenient to choose a smaller reference 
length scale when one is interested in short distance properties of the theory. We 
therefore introduce also 

rc = r(0.65) (2.2) 
and would like to know its relation to tq. 

2.2 The ratio Tc/tq 

The evaluation of r(c) from eq. ( ^.Ij ) begins with the extraction of the static quark 
potential V{r) between quarks separated by a distance r along a lattice axis. We 
use Wilson loop correlation functions from as well as new simulations for 
smaller values of the lattice spacing (/5 = 6.57, 6.69, 6.81, 6.92) keeping the size 
of the L'^-torus at L 3.3ro, where finite size effects in V{r) are known to be 
small for r < rg. For the smallest lattice spacing this required a 64^ lattice. We 
essentially followed the procedure of [|T6|. Some details are given in App. 0. 
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Figure 1: The ratio rc/ro for 5.95 < (3 < 6.57 in the standard Wilson action (circles) including 
the continuum extrapolation (solid point). This is an analysis of data in |]l6| . 



The force at finite lattice spacing is then defined as 

F(ri) = [V{r) - V{r - a)] /a, (2.3) 

with the distance ri chosen such that the force on the lattice has no deviations 
from the force in the continuum when evaluated at tree level. Lattice artifacts 
are then suppressed by a power of a; F{ri) is a tree-level improved observable. In 
the SU(2) theory it has been observed that the remaining lattice artifacts in the 
force are surprisingly small and here we shall again find no evidence for them 
as long as r > 2a. 

Explicitly we have 

{Anr^)-^ = [G{r, 0, 0) - G{r - a, 0, 0)] /a, (2.4) 

where G{r) is the (scalar) lattice propagator in 3 dimensions defined in eq. ([B.!]). 



With the methods of |T8[ it can be constructed directly in coordinate space. To 
our knowledge these methods have not been applied in 3-d and we had to find 
the necessary invariant of the recursion relation eq. ( p.3|) . Formulae are listed in 
App-l 

Solving eq. ( p.l|) requires furthermore an interpolation of F{ri). This can 
be done with small systematic errors [Q, which nevertheless dominate over the 
statistical errors at small r (see App. ^for details). 

Figure |1| shows the ratio Tc/to for several lattice spacings. No dependence on 
the resolution is seen within the errors of below 1%. A continuum extrapolation 
gives 

rjro = 0.5133(24) (2.5) 
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Figure 2: Interpolation of ro/a. 



and we note that it is safe to use eq. ( p.5D also at finite lattice spacings starting 
around P = 6.4. 

2.3 Parameterization of ro/a 

The direct determination of ro/a for 5.7 < (3 < 6.4 [|I6| Q and our new computa- 
tions of Tc/ a in the range 6.57 < {3 < 6.92 may be combined with rd ro=0. 5133(24) 
to obtain an interpolating formula giving ro/a in the whole range 5.7 < (3 < 6.92. 



Following |jT6| we interpolate ln(a/ro) through a polynomial in (3 and find that 



ln(a/ro) = -1.6804 - 1.7331(/3 - 6) + 0.7849(/3 - 6)^ - 0.4428(/5 - 6)^ 

for h.l <I3< 6.92, (2.6) 

is an excellent approximation (Fig. Q) to the MC results. This formula covers a 



larger range of [3 than the one given in |jl6|, but its precision in the low 13 range 



is somewhat worse. The accuracy of ro/a in eq. ( p.6| ) is about 0.5% at low (3 
decreasing to 1% aX (3 = 6.92. 

^ In fl^ the simulation for /3 = 6.57 was performed with a single smearing level, yielding less 
control over excited state contaminations. As a result the error on ro/a seems to be somewhat 
underestimated. Although the statistical significance of this small effect is not clear, we decided 
to use our new data for (3 = 6.57 instead. 
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Table 1: Results for r^/a ||T^ and Vc/a. 



/3 


ro/a 


/5 




5.7 


2.922(9) 


6.57 


6.25(4) 


5.8 


3.673(5) 


6.69 


7.29(5) 


5.95 


4.898(12) 


6.81 


8.49(5) 


6.07 


6.033(17) 


6.92 


9.82(6) 


6.2 


7.380(26) 






6.4 


9.74(5) 







3 Continuum force and potential 



We now determine the continuum force and potential by extrapolation of the 
MC-results at finite values of the lattice spacing to the continuum. It is expected 
that close to the continuum limit the dominant discretization error is quadratic 
in the lattice spacing. A clean argument why this is so has never been given in 
the literature. We shall fill this gap and give one in the next section. 

Figure § shows the dependence of the force in units of on the resolution 
for some selected values of the separation r. Fitting 



rlF{r) = rlF{r) _ 1 + s x (a/rj 



a=0 



(3.1) 



for fixed r/rc, the slope s is statistically not significant throughout our range of r 
and j3. Our statistical precision allows to quote 



|s| < 1 for r/a > 2 , 0.4 < r/rc < 1 • 



(3.2) 



Of course this 1-cr bound is valid only when all details are as discussed above, in 
particular eqs. ( ^.3| , ^.4D are used to define the force at finite a. If the naive form 
ri = r — I is employed instead, the corresponding slopes s become rather large, 
as can be seen in the figure. 
A similar statement. 



tIF{t) 



rlF{r) 



a=0 



s X {a/roY 



\s\ < 1.2 for r/a > 2 , 0.5 < r/ro < 1.5 , 



(3.3) 
(3.4) 



can be made for larger r. 

The continuum force is plotted in Fig. |^ using eq. (|2.5| ) to combine the two 
regimes of r. Some data at finite /3 are included in the figure. In these cases we 
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Figure 3: Continuum extrapolation of r'^F{xrc), for x — 0.4, 0.5, 0.9 from top to bottom and 
oi rQF{xro)j for x — 0.5, 0.6, 1.5 from top to bottom. The data are from our new computations 
and from jl^. Filled circles correspond to the naive value ri = r — | instead of eq. (2.4). 



used our "bounds" on s to estimate that the discretization errors are smaller than 
the statistical ones. 

For large values of r, the force is expected to be given by a constant, the string 
tension, plus a first universal [|1^ correction. Assuming this description to 
be valid already at r = rg yields the parameter free bosonic string model, 



F(r) =cr + 



12r2 



arl = 1.65-7r/12, 



(3.5) 



which is in excellent agreement with our results for r > 0.8ro. Note that in the 
same region of r, excited potentials do not at all follow the expectations from 
an effective bosonic string theory [Q. This suggests that the agreement with 
eq. ( p.5|) is rather accidental. In any case one would expect corrections to this 
formula to be negligible only for much larger r. Nevertheless eq. ( p.5|) is a very 
good effective description of F{r) for 0.8ro < r < l.Gro- 

At short distances the force may be obtained by an integration of the pertur- 
bative renormalization group. 



F{r) = CFgl^{r)/{4nr'), = 4/3 
d _ 



bo 



167r2 ' "1 



(3.6) 
(3.7) 



u=0 

102 
(167r2)2 



&2 = (-3470 + 2519f - 99^ + 726C(3); 
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Figure 4: The force in the continuum limit and for finite resolution, where the discretization 
errors are estimated to be smaller than the statistical errors. The full line is the perturbative 
prediction with Aj^ro — 0.602. The dashed curve corresponds to the bosonic string model 
normalized by r'gi^(ro) — 1.65. 



Here the 3-loop coefficient, 62, could be extracted from |7|j6|j2l|]. Inserting |l|] 

Aj^ro = 0.602(48) (3.^ 



as well as the known relation between and Aqq to fix the integration 
constant, we have a parameter free perturbative prediction. It agrees with the 
non-perturbative results up to r «i O.Sro. Indeed, inserting Aj^tq at the upper 
end of the error bar of eq. ( p.8|) into the perturbative formula, very close agreement 
with the data points is seen in this range of r. Again this agreement at quite 
large r appears somewhat accidental as the perturbative prediction itself is only 
stable at smaller distances. Certainly there is no need for large non-perturbative 
terms at such distances as it was concluded earlier on the basis of an exploratory 



investigation p2[. We shall discuss the comparison with perturbation theory in 
more detail in a separate publication. 

The potential contains the same physical information as the force but statis- 
tical and systematic errors in the lattice determination are different. To eliminate 
the self energy contribution we consider Vi(r) — li(rc), where we apply a tree-level 



7 



1.8 



1.7 



1.6 
0.21 



> 0.2 

X 0.19 
X -0.02 

-0.03 

-0.04 



I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' M 

1 -0.067 

-0.0675 

-0.068 

1- 1 I I I I I I I I I I I I I I I I I I I I I I J -0.0685 





I I I I I I I I I I I I I I I I H 
0.02 0.04 0.06 0.08 




0.01 0.02 0.03 0.04 

(a/rj^ 



0.01 0.02 0.03 

(a/rj^ 



Figure 5: Continuum extrapolation of Vi{r), for r/ra = 1.5, 0.6, 0.5 on tlie left hand side and 
for r/rc — 0.9, 0.4, 0.3 on the right hand side. 



improvement as for the force, 

Vi{di) = V{d) , (Arrdi)-^ = G{d, 0,0). (3.9) 

The continuum extrapolation is performed exactly as for the force. In this case 
slopes, s, defined as above are statistically significant for both our smallest dis- 
tances (r ~ 0.3rc) and the largest one (r ^ 1.5ro). As a consequence, at small r 
the combination [Vi{r) — Vi{rc)]rc has about 1% discretization errors ai jS = 6.92 
while in the large distance region these errors go up to 1% in [Vi(r) — Vi(rc)]ro 
at P = 6.4. These are the /5-values corresponding to the smallest lattice spacings 
available in the two different regions. 

The continuum potential is plotted in Fig. ^, with some data at finite f3. 
In this figure the statistical as well as the discretization errors are below the 
size of the symbols. In the whole range of r the non-perturbative results are 
described by the model eq. (|3.5| ) within about 1% accuracy. For short distances, 
r < 0.3ro, the perturbative prediction V{r) = V{0.3rc) + Jo^^rc ^V^iv) (with the 
perturbative expansion for F as discussed above) is quite accurate. For future 
reference continuum force and potential are listed in Table ^ 
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Figure 6: The static potential. The dashed Une represents the bosonic string model and the 
solid line the prediction of perturbation theory as detailed in the text. 
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Table 2: Potential and force after continuum extrapolation. 



r/fc 




rQ{vi[r) — vi[rc)) 




r^t [r) 




0.6 
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n A 
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u.u 




u. iyoy ( u J 


0.5 


5.363(78) 


-0.8926(74) 


0.7 


1.951(14) 


0.4051(89) 


0.6 


4.244(53) 


-0.6475(54) 


0.8 


1.812(11) 


0.5930(99) 


0.7 


3.538(48) 


-0.4494(35) 


0.9 


1.722(11) 


0.769(11) 


0.8 


3.060(38) 


-0.27950(15) 


1.1 


1.592(10) 


1.101(11) 


0.9 


2.713(30) 


-0.13259(67) 


1.2 


1.559(13) 


1.258(12) 








1.3 


1.537(18) 


1.413(13) 








1.4 


1.537(24) 


1.567(14) 








1.5 


1.507(27) 


1.722(24) 



4 Lattice artefacts 

The standard framework for the discussion of lattice artefacts is Symanzik's effec- 
tive theory, which is expected to give the asymptotic expansion of suitable lattice 
observables in integer powers of the lattice spacing (up to logarithmic modifica- 
tions) in an asymptotically free theory In the 2-dimensional 0(3) a model 



it predicts this expansion to start at order a? but unexpectedly numerical results 
are described better by a dominant linear term in a for a range of lattice spacings 
I13. While there is no evident contradiction with the result of an analysis a la 



Symanzik which is supposed to describe the asymptotic behavior, the standard 
picture should be tested in 4-d gauge theories and QCD as much as possible. 

Below we shall consider some observables and study their a-dependence. Be- 
fore coming to these numerical examples, we want to give a brief but thorough 
argument that the leading artefacts in potential differences are expected to be 
O(a^). This is not obvious since the potential is usually defined in terms of a 
Wilson loop, which does not fall into the category of correlation functions of local 
fields discussed by Symanzik. 

4.1 Why are the leading lattice spacing errors to F{r) quadratic in a? 

In a short version, our answer to this question is that the heavy quark effective 
theory |^ formulated on the lattice is 0(a)-improved without adding any addi- 
tional operators to the Lagrangian The potential is an energy of the effective 
theory and thus has no linear a-effects. Of course, for this statement to be mean- 
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ingful, the self energy has to be ehminated by considering potential differences or 
the force. We want to also mention that our argument is based on the assumption 
that the effective theory is renormalized as usual by the addition of local opera- 
tors. Explicit perturbative computations support this assumption but it has not 
been proven so far. 

We find it easiest to discuss the a-effects starting from a correlation function 
in the effective theory with Schrodinger functional boundary conditions p7| , |28| , p6| . 
A correlation function is chosen which describes a quark-antiquark pair separated 
by X in space and propagating in Euclidean time from xq = to Xq = T. Boundary 
conditions are taken exactly as specified in the quoted papers (C = C" = 0). In 
the notation of p6|, the correlation function is given by 



A;hh(x,T) = (Ch(0)75Ch(x) Ch'(x)75Ch'(0)) . (4.1) 
In the standard formulation of the lattice theory, a positive transfer matrix exists 



29 1 and therefore fchh has the exact representation 



fchh(x, T) = \Y: B^e'^-^ Cn(x)e-^" (4.2) 

with positive B^.Cn^'^)- The energies are the (negative) logarithms of the 
transfer matrix in the vacuum sector and ^(x) are the energies in the proper 
charge sector (see e.g. for a discussion of the transfer matrix in the presence 



of static charges). Adopting the convention Eq = Q (vanishing vacuum energy), 
the potential is given by V{r) = Vo(r, 0, 0). 

Renormalization and improvement of /chh(x, T) follows directly from the dis- 
cussion of [2C]: the renormalized and improved correlation function, 



fchh,R(x, T) = e-2^™^ A;hh(x, T) , (4.3) 



satisfies 



x,T)= VR(x,T)L^o + 0(a2) (4.4) 



and the renormalization constants 6m and Zh may be chosen to depend on nothing 
but the bare coupling. From eq. (|4.4| ) we may conclude 



Vin) - Vir,) = [Vin) - Vir,)l^, + Oia') , (4.5) 



^ It is assumed here that the improvement coefficient ct is chosen properly to remove 0(a) 
effects specific to the Schrodinger functional pl[. When full QCD with quarks is considered, 



one furthermore has to account for Csw,ct [p2[- In massive QCD, the term proportional to &h 



has to be included in eq. (4.3) but plays no role in the following. 
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since 



and 



V{r) = Voir, 0, 0) = a'^ lim In W(x,T a) 



Our result eq. (O) is valid for the pure gauge theory as well as for the 0(a)- 



improved formulation of QCD discussed in |33,|32|. 



4.2 Lattice spacing effects in force and potential. 

Examples for a-effects are shown in Fig. |^ and Fig. ^. The following observations 
are relevant. Using eqs. (|2.3|J2.4| ), i.e. a tree-level improved definition of the force, 
lattice spacing effects are below the level of our small statistical errors when one 
restricts oneself to r > 2a and the lattice spacings considered here. By contrast, 
with a naive midpoint rule for the force, the a-effects are quite sizeable. The 
difference is by construction a pure power series in a, starting with a^. For r > 2a 
this series is well approximated by the leading term. Thus both data sets in Fig. | 
contain the same essential information: full compatibility with Symanzik's theory 
of discretization errors. The case of finite potential differences. Fig. differs only 
slightly. Here, small a-effects could also be observed for the tree-level improved 
definition. Again - but not shown in the figure - the standard definition without 
tree-level improvement shows quite large a-effects. This has been known for a 



long time p4| , p!0| . In the past, usually a correction term, relying on a fit to the 
potential, has been applied to remove the dominant a-effects [^. As shown in 
PI , and again here, this procedure may be replaced by a simple and theoretically 
sound definition. 

4.3 Continuum extrapolation of Lmax/'^o 

In a calculation of the A-parameter of the pure gauge theory was presented. 



The result may be split in the following way, 

Aj^ro = (2.0487AL„,,J/(L^ax/ro) • (4.^ 



Here L^ax is the distance where the Schrodinger functional coupling |^ has the 
value g^{L^a.x) = 3.48 and the combination AL^ax = 0.211(16) in the first paren- 
thesis is obtained from the non-perturbative evolution of the Schrodinger func- 



tional coupling. The ratio -Lmax/^'o "was computed in [16 



= 0.718(16). (4.9) 

a=0 
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Figure 7: Fits to eq. ( 4.10 ). The fit represented by the full lines assumes p"^ '""^^ — 0, while 
for the dotted lines all coefficients are left as free parameters. 



Now we know tq closer to the continuum and can further investigate the contin- 
uum extrapolation which lead to eq. ( [4 .91) . Details of the calculation, which are 
described in [|16|, are not repeated here. We only recall one important point. For 
the Schrodinger functional of the pure gauge theory, Symanzik's analysis implies 
that lattice artifacts linear in the lattice spacing exist in general. They should 
vanish when the improvement coefficient Ct (specific to the Schrodinger functional 
) is chosen properly. Since Ct is know to 2-loop precision, it is interesting to com- 
pare results for L^^ax/ro obtained with Ct at 2-loop accuracy to those using Ct at 
1-loop. 

We obtained new data points by adding an entry L/a = 14 , P = 6.926(5) to 
Table 3 of [|16]| and using Vc/a and Vc/ro from Sect. ||. The resulting ratios are 
shown in Fig. ^ together with fits of the form 



To 



ro 



+ Pi — +P2 — , 

a=o '^o 



1-loop, 2-loop . 



(2-loop) 



(4.10) 



0. Our 



If 2-loop approximation for Ct is sufficient, we expect to have pi 
results are compatible with this assumption (full lines, full square). On the other 
hand, when we leave p"i as a free parameter, the best fit contains a noticeable 
linear term and the continuum point (i^max/''^o)a=o differs significantly (dotted line 
and open square). Despite the precision of better than 1% and the considerable 
range in the lattice spacing, our data are not enough to decide between the two 

We thank Jochen Heitger for performing the necessary simulation of the Schrodinger func- 
tional . 
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types of fits, which both have a good x^- As continuum hmit, we quote 



To 



= 0.738(16), (4.11) 

a=0 



covering the range allowed by both fits. Due to the uncertainty in the proper 
extrapolation formula, the error is not reduced compared to eq. ( |4.9| ). The small 
change in the central value hardly affects the value for the A-parameter. We now 
have Aj;js = 0.586(48)/ro compared to A^g = 0.602(48)/ro of 0. 

Our inability to improve the precision of the continuum limit eq. ( [4.11D does 



not signal a problem with the standard theory of lattice artefacts. Rather we 
do have a difficult situation here where both linear and quadratic lattice spacing 
effects are expected and indeed both appear to be significant. 



5 Discussion 

In the pure gauge theory it was possible to compute the static potential over a 
large range of distances by considering two sets of lattices. Very small lattice 
spacings were used in the short distance range and larger ones for r > 0.2fm. In a 
region of overlap the two sets of computations could be matched (in the continuum 
limit). We find that the continuum force Fig. ^ is a nice demonstration of the 
potential of lattice simulations for precision physics. The force is in satisfactory 
agreement with perturbative predictions. A non-trivial point in this comparison 
is that the perturbative expression does not involve any free parameter since the 
A-parameter in units of tq is known. Nevertheless there is a question of the 
truncation of the perturbative series. Here we truncated the /3-function at 3-loop 
order. Other possibilities and the stability of the perturbative prediction will be 
discussed in a separate publication. 

We have also investigated lattice spacing effects in some detail. First of all we 
have considered the question whether close to the continuum the leading lattice 
artefacts are quadratic in a as commonly has been assumed. A positive answer to 
this question could be given by formulating it in terms of the heavy quark effective 
theory for which Symanzik's discussion of a-effects is expected to apply. Also our 
numerical results are in full agreement with O(a^) lattice artifacts in potential 
differences. 

A numerically difficult situation is met when linear and quadratic a-effects 
compete. This is the case in the ratio L^a^/r^ considered in Sect. ^]3[ Here the 
removal of linear a effects is incomplete since the improvement coefficient Ct is 
known only perturbatively. Despite accurate numbers at finite lattice spacing our 
continuum estimate eq. ( [4. 11] ) has a relatively large uncertainty. In order to avoid 
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such a situation in simulations of full QCD, non-perturbative 0(a)-improvement 
p6| should be applied or one should use formulations with exact chiral symme- 
try where linear terms are absent without any necessity of tuning improvement 
coefficients |p7| . 

We finally point out that at our level of precision the observed lattice spacing 
effects are in complete agreement with the standard theory of Symanzik. Our 
precision does however not quite reach the one achieved in 2 dimensions, where 
indications for deviations from the standard picture have been found . 
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A Computation of potential and force 

We computed Wilson loop correlation matrices using smearing exactly as in |jl6| 



Three different smearing levels n;, / = 0, 1, 2 were used, constructing a 3 x 3 matrix 
for each value of the loop size r, t with a < r < r^ax and a <t < tmax- Our choices 
for the different parameters are listed in Table 0. The numerical values for rii, 
defined as in |jT6[, are in general smaller than in that reference since we want to 



compute the potential for shorter distances. For the same reason only 40 to 75 
measurements were sufficient to get satisfactory precision. With relatively modest 
computational effort we could extend the calculations to a ^ 0.025fm where a 64^ 
lattice had to be simulated. Several iterations of a hybrid overrelaxation algorithm 
with Nor overrelaxation sweeps per heatbath sweep were performed to separate 
the field configurations of our MC sample. 

The potential V{r) was extracted from the correlation matrices following 
the procedure of []TB|. As in this reference, we found that for t > excited state 
contaminations are negligible when exp[— to(^i(^) < 0.3 is satisfied. Since 

the gap Vi{r) — Vo(r) grows towards small r from 3/ro at r = Tq to ~ 7/ro at 
r = O.lro, the short distance region is easiest in this respect. All errors quoted in 
this paper were computed by jacknife binning. Force and potential are listed in 
Tables 4-7. 
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Table 3: Simulation parameters. 



L/a 


/9 


no 


rii 


n2 


Nor 


'"max/ Oi 


^max/ 


40 


6.57 


0, 


77, 


153 


18 


7 


11 


48 


6.69 


0, 


53, 


106 


22 


10 


12 


56 


6.81 


0, 


72, 


144 


25 


12 


14 


64 


6.92 


0, 


94, 


188 


29 


12 
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Interpolation. 

Solving eq. ( p.l|) and evaluating force and potential at distances r = xr^ for 
given X, requires their interpolation. While this can in principle be done in many 
ways, it is advantageous to use physics motivated interpolation formulae. 

For the force we followed choosing the interpolation function 

F{r) = f, + f2r-', (A.l) 

between the two neighboring points. The systematic error arising from the interpo- 
lation was estimated adding a term /sr"^ and taking a third point. The difference 
between the two interpolations was added (linearly) to the statistical uncertainty. 
We observed that at least for r^0.4ro the interpolation error is smaller than the 
statistical uncertainty. For small distances the systematic error increases and can 
be of the order of the statistical one or even bigger. 

The potential is interpolated by the corresponding ansatz 

Vi{r) = vi + V2r + vsr'^ . (A.2) 

Here two points are chosen such that the desired value of r is in between. For 
the choice of a neighboring third point one has two possibilities, leading to two 
results. Their difference was taken as the interpolation error. Also in this case 
the systematic errors are larger than the statistical ones at short distances while 
at large distances the situation is reversed. 

B The evaluation of the 3-dimensional lattice propagator 
in coordinate space 

An efficient method to calculate the lattice propagator in coordinate space is 



proposed in |38[. It is based on a recursion relation which allows to express the 
propagator as linear function of its values near the origin. In this paper, the 4- 
dimensional case is discussed; it is straitghforward to apply the same method in 
3 dimensions. We use lattice units, a = 1, in this appendix. 
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Starting from the Laplace equation for the Green function, 

- AGfx) = I '^f = ° , (B.l) 
otherwise 

with 

A/(x) = Y: [/(x + j) - 2/(x) + /(x - ])] , (B.2) 
and using other general features of the lattice propagator, one obtains (j = 1, 2, 3) 

3 

Glx + j) = G(x - j) + 2^ J2[G{^) - - ^)] , (B.3) 

P ^=l 

forp = E?=ia;, ^0. 

Because of isotropy, eq. ( [B.3| ) can be restricted to the points x with xi > 
X2 ^ X3 > 0. In this region eq. ( [B.3| ) can be used as recursion relation to express 
G(x) as a linear combination of G(0, 0, 0), 0, 0), 1, 0), 1, 1). These 
four values at the corners of the unit cube are not independent: from eq. ( |B.1| ) at 
x= 

G(0,0,0)-G(1,0,0) = | (B.4) 

follows directly. Another relation can be deduced in the following way: first one 
observes that eq. (|B.3|) becomes one-dimensional along the lattice axes. Defining 

gi{n) = G{n, 0, 0), g^in) = G{n, 1, 0), gs{n) = G{n, 1, 1) (B.5) 

and using the lattice symmetries, we find 

gi{n + l) = 6(7i(n)-%2H-(7i(n-l), (B.6) 
92{n + l) = ^A^92{n) - g^{n) ~ g,{n)] ~ ^^g,{n - 1) , (B.7) 
93{n + l) = ;g2[3^3(n)-2^2(n)]-^^3(n-l). (B.8) 

Next one notices that 

k{n) = {n - l)gi{n) + 2ng2{n) + {n + l)gs{n) - ngi{n - 1) - (B.9) 
2(n - l)(?2(n - 1) - (n - 2)(73(n - 1) 

is an invariant of the recursion relation, that is k{n + 1) = k{n) for n > 1. The 
value of k is worked out in the limit n 00, where 

9A^) = ^ + 0{l/n'), for j = 1,2,3. 
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This consideration yields k{n) = for n > 1. Setting n = 1 in eq. ( [B.9| ) one 
obtains 



3G(1, 1, 0) + 2^(1, 1, 1) - ^(0, 0, 0) = 0. (B.IO) 

Thus from the four initial values two can be eliminated and the propagator is 
obtained in the form 



G(x) = ri(x)G(0, 0, 0) + r2(x)G(l, 1, 0) + r3(x). (B.ll) 

The coefficients ri,r2,r3 are rational numbers and can be evaluated recursively 
from eq. ( |B.3| ). 

From the numerical point of view, we are then left with the task to accurately 
compute G{0, 0, 0) and 1, 0). This can be done with the procedure discussed 
m givme; 



G(0,0,0) = 0.2527310098586630030260020266135701299..., 
G(1,1,0) = 0.0551914336877373170165449460300639378... 
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Table 4: Force and potential in the short distance region. 



B 


1 1/ Uj 


u r yi ij 




n V,(di) 








1 855 


n 4578Q8f71 1 




9 977 




9 RRQ 






3.312 


0.033391(97) 


3.922 


0.54795(21) 




4.359 


0.023752(94) 


4.942 


0.57170(28) 




5.393 


0.01904(12) 


5.954 


0.59074(35) 




6.414 


0.01629(13) 


6.962 


0.60703(41) 


R RQ 








n 4'^Q18f41 




9 977 






n 4Q1 AR7(70\ 




SI 9 




S Q99 








n 091 n^zLi'7^'1 


A Q49 






O.Oc/O 


n niR4'^Q('79l 


5 Q54 






414 


U.UXO / ZiOIOZi 1 


R QR9 


^79R4r9R'i 






ni9n'^RfQR'l 

U.UXZiUOUl t7U J 


7 QR7 


U.U04:0/ I OU J 




8 438 


n mn8fi9('82'i 


8 971 


n 5Q574('41 1 




9.445 


0.010123(97) 


9.974 


0.60587(46) 


6.81 






1.855 


0.422617(22) 




2.277 


0.048833(26) 


2.889 


0.471411(47) 




3.312 


0.027650(32) 


3.922 


0.499062(66) 




4.359 


0.018860(31) 


4.942 


0.517921(90) 




5.393 


0.014471(32) 


5.954 


0.53239(11) 




6.414 


0.011870(40) 


6.962 


0.54426(13) 




7.428 


0.010163(54) 


7.967 


0.55437(18) 




8.438 


0.009072(56) 


8.971 


0.56344(21) 




9.445 


0.008267(52) 


9.974 


0.57171(24) 




10.451 


0.007701(58) 


10.977 


0.57941(27) 




11.455 


0.007232(62) 


11.979 


0.58664(30) 
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Table 5: Force and potential in the short distance region. 





ri/a 


F{ri) 


di/a 


a Vi{di) 


6.92 






1.855 


0.408642(19) 




2.277 


0.046081(19) 


2.889 


0.454723(34) 




3.312 


0.025696(24) 


3.922 


0.480418(47) 




4.359 


0.017266(29) 


4.942 


0.497662(75) 




5.393 


0.012969(33) 


5.954 


0.510631(88) 




6.414 


0.010412(39) 


6.962 


0.52104(11) 




7.428 


0.008855(37) 


7.967 


0.52990(13) 




8.438 


0.007755(43) 


8.971 


0.53765(16) 




9.445 


0.006974(51) 


9.974 


0.54463(19) 




10.451 


0.006402(53) 


10.977 


0.55103(23) 




11.455 


0.006022(53) 


11.979 


0.55705(26) 
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Table 6: Force and potential from the data of Guagnelli et al. 



(3 


n/a 


F{ri) 


di/a 


a Vi{di) 


5.95 






1.855 


0.61794(16) 




2.277 


0.11212(20) 


2.889 


0.73006(33) 




3.312 


0.08319(27) 


3.922 


0.81325(55) 




4.359 


0.07164(36) 


4.942 


0.88489(84) 




5.393 


0.06613(48) 


5.954 


0.9510(12) 




n ATA 

6.414 


0.06296(56) 


6.962 


1.0140(17) 




7.428 


0.0606(12) 


7.967 


1.0728(30) 


6.07 






1.855 


0.571729(97) 




2.277 


0.09211(11) 


2.889 


0.66384(19) 




3.312 


0.06427(13) 


3.922 


0.72811(30) 




4.359 


0.05301(26) 


4.942 


0.78116(55) 




5.393 


0.04771(28) 


5.954 


0.82887(72) 




6.414 


0.04468(27) 


6.962 


0.87355(89) 




7.428 


0.04262(38) 


7.967 


0.9162(11) 




8.438 


0.04215(45) 


8.971 


0.9583(13) 




9.445 


0.04087(73) 


9.974 


0.9992(17) 


6.2 






1.855 


0.533457(82) 




2.277 


0.07804(11) 


2.889 


0.61145(17) 




3.312 


0.05135(14) 


3.922 


0.66279(28) 




4.359 


0.04054(13) 


4.942 


0.70333(38) 




5.393 


0.03511(20) 


5.954 


0.73844(54) 




6.414 


0.03238(20) 


6.962 


0.77082(69) 




7.428 


0.03018(25) 


7.967 


0.80100(85) 




8.438 


0.02884(26) 


8.971 


0.8298(10) 




9.445 


0.02813(27) 


9.974 


0.8580(12) 




10.451 


0.02766(30) 


10.977 


0.8856(14) 




11.455 


0.02752(34) 


11.979 


0.9131(16) 
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Table 7: Force and potential from the data of Guagnelli et al. ||T^ 





ri/a 


F{ri) 


di/a 


a Vi{di) 


6.4 






1.855 


0.488379(48) 




2.277 


0.064318(51) 


2.889 


0.552697(81) 




3.312 


0.039580(70) 


3.922 


0.59228(13) 




4.359 


0.029360(79) 


4.942 


0.62164(19) 




5.393 


0.024367(84) 


5.954 


0.64600(24) 




6.414 


0.02145(12) 


6.962 


0.66744(40) 




7.428 


0.01939(16) 


7.967 


0.68683(50) 




8.438 


0.01819(18) 


8.971 


0.70502(66) 




9.445 


0.01757(17) 


9.974 


0.72258(79) 




10.451 


0.01677(17) 


10.977 


0.73936(84) 




11.455 


0.01651(15) 


11.979 


0.75586(94) 




12.459 


0.01609(17) 


12.980 


0.7720(11) 




13.462 


0.01616(29) 


13.982 


0.7881(11) 




14.465 


0.01564(18) 


14.983 


0.8038(12) 




15.467 


0.01513(39) 


15.984 


0.8189(14) 
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